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STATISTICAL  CORRELATION  DISTANCE  EOR 
PLATFORM-TO-PLATFORM  RADAR  INTEGRATION 


INTRODUCTION 

In  search-radar  tracking  systems  the  detections  must  be  associated  or  correlated  with  the 
tracks.  One  way  of  correlating  the  detections  with  tracks  is  to  compute  a statistical  distance 
between  each  track  location  and  measurement  and  choose  the  detection  to  go  with  the  track 
which  has  the  minimum  distance  function.  Usually  those  detection  and  track  combinations  are 
limited  to  a small  region  of  space  such  that  the  distance  would  be  below  some  threshold,  thus 
reducing  the  number  of  calculations  required  dramatically.  This  type  of  correlation  has  been 
used  with  most  search-radar  tracking  systems  to  date. 

This  report  will  first  rederive  the  basic  statistical  distance  used  in  most  collocated  radar 
tracking  systems.  The  derivation  will  be  based  on  the  standard  hypothesis-test  procedure.  Next 
the  test  will  be  complicated  somewhat,  and  a new  statistical  measure  will  be  obtained  which 
would  probably  be  of  more  use  in  integrating  radar  data  from  noncollocated  sites. 

STANDARD  STATISTIC  AL  DISTANCE 


We  begin  formulating  the  standard  hypothesis  test.  Given  a detection  and  measured  posi- 
tion Xm,  find  which  hypothesis  Hr  i = 1 M , has  the  highest  probability  of  being  true.  This 

is  given  by 


P(Hj/Xm  ) > P(HJXm  ), 


(1) 


where  H,  is  the  hypothesis  that  the  detection  belongs  to  /th  track.  This  should  be  read  that  the 
hypothesis  Hj  is  chosen  if  the  probability  of  being  true  given  the  measurement  is  greater 
than  the  probability  of  //,  being  true  given  the  measurement  for  all  values  of  /.  By  use  of 
Bayes  rule  the  hypothesis  test  can  be  rewritten 


P(Xm/Hj  )P(Hj ) P(XJH,)PW,) 
P(Xm)  ^ P(XJ 


(2) 


The  P(Xm ) is  irrelevant  and  can  be  removed.  The  measured  position  can  be  written  as 

Xnl  = Xr  + Nm,  (3) 

where  XT  is  the  target  true  position  and  Mm  is  the  measurement  noise.  By  change  of  variable 
using  (3),  the  hypothesis  test  (2)  can  be  rewritten  as 


PK  (Xm  - Xt  /H,  )P(Hi  ) 3s  PNm  (Xm  - XT  /H,  )P(H, ). 

since  the  noise  and  the  hypothesis  are  independent,  (4)  becomes 

/><//,) 


P.  < xm  -xT)>  Ps  ( Xm  - XT  ) 


P(H, ) 


(4) 


(5) 
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B II  CAN  TRIiLL 


If  the  noise  is  Gaussian  or  near  Gaussian,  the  hypothesis  test  (5)  is  equivalent  to 

[U„,  -XT)  - (Xm  - XT)V  RN-J  uxm  — XT)  - (Xm  -XT)\ 

< K*w  - XT  ) - (Xm  - XT  )\  Rs-'  ixm  - X,  ) - (X„,  - Xj  ) ] + C,  (6) 

w here  R ^ l is  the  covariance  matrix  of  noise,  C — In  PiHJPiH, ),  and  the  prime  denotes  the 
transpose.  If  the  measurement  is  unbiased  (Xm  = XT  ) and  the  hypotheses  have  equal  proba- 
bility of  occurrences  (P(HI ) = P(Hf ) ),  then  (6)  reduces  to 

1 iX'»  ~ Xr,  > ) K yj  t - XTl  > 1 < I ixm  -Xr))  I (Xm  -XT)l  (7) 

For  "good"  collocated  radars  the  measurements  can  be  made  unbiased  in  most  cases.  (This  is 
not  true  in  gener  I unless  considerable  effort  at  the  radar  level  is  made  and  bias-removal  algo- 
rithms are  used  between  different  collocated  radars.)  With  the  assumption  of  no  bias,  for  collo- 
cated radars  the  quantities  are  defined 

v 1 v 1 o -i  _[ l/,rR  o I 


y m 

A"'  ~ a:„ 


XT  ~ ah 


, and  R„ 


The  noise  is  uncorrelated  in  the  range  and  azimuth  measurement,  and  the  hypothesis  test  be- 


rr  ) ‘ 


( a:,„  — a:T  ) 2 


(r„,  - rT  ) 2 


{a:m  ~ a:T,  )' 


It  is  assumed  that  the  ir K and  <ra:  are  estimated  when  the  measurement  is  taken,  and  rm  and 
a:m  are  the  measured  quantities.  However,  the  true  position  rT  and  a:T  are  not  known.  To 
complete  the  test,  rT  and  azT  must  also  be  estimated.  The  best  estimate  of  the  target's  posi- 
tion comes  from  the  tracking  filter,  which  may  be  the  Kalman  filter  or  some  variation  of  it.  In 
any  case  the  best  estimate  of  the  target's  true  position  is  the  predicted  position  of  the  filter 

x„  = XT  + Nn  = _P_  1.  (10) 


Using  the  best  estimate  of  the  target’s  position  XP  , the  best  estimate  of  the  measurement  error 
it R and  <ra:.  and  the  measured  position  Xnn  the  hypothesis  test  becomes 


{rm  ~ V > 2 


(az,„  — azP  ) ‘ 


( az„ , — azi 


The  statistical  distance  D,  is  defined  as 


(r„,  ~ rP  )2 


(azm  - azp  ) 2 


This  is  the  usual  statistical  distance  used  in  minimum-distance  correlators  for  tracking  with  col- 
located radars.  The  statistical  distance  is  also  sometimes  used  for  determining  if  a target  is  turn- 
ing or  not.  For  example,  if  the  predicted  position  is  biased  from  the  true  position,  the  distance 
D becomes  large  and  the  bandwidth  of  filter  should  be  opened  up  if  Ihe  association  was 
correct.  At  this  point  we  will  consider  a somewhat  more  complicated  correlator  based  on 
knowing  both  the  measured  and  predicted  position. 
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NEW  STATISTICAL  DISTANCE 

The  formulation  of  the  new  statistical  distance  parallels  the  previous  development.  The 
hypothesis  is 

P(HtIXm.Xf)  > P(HJXnr  Xp  ),  (13) 

in  which  //,  is  the  hypothesis  that  the  measurement  Xm  belongs  to  the  /th  track, 

xm  = x,  + Nm-  (14) 

and 

Xp  -X,  + Np.  (15) 

where  Xm  is  the  measured  position,  X , is  the  target’s  true  position,  Nm  is  the  measurement 
noise,  and  Np  is  the  predicted  position  noise.  By  use  of  Bayes  rule  the  hypothesis  test  (13)  can 
be  written  as 

P(Xm.Xp/HJ)P(HJ)  > P(Xm.Xp/H)P(Hl)  (]6) 


P(Xm,Xp) 


P(Xm.Xp) 


The  P(Xm,  Xp)  is  irrelevant  to  the  test  and  can  be  ignored.  Equations  (14)  and  (15)  can  be 
used  to  change  the  variable  in  (16),  yielding 

PNm.N,  1 <Xm  - Xt/  ),  (Xp  - X,'  )/HJ  )P(Hj ) 1 

3s  PNmiN'  HXm  - X,  ),  (Xp  - X,  )/H,  )P(H, ) ].  (17) 

The  noise  is  independent  of  the  hypothesis,  and  the  a-priori  probabilities  are  assumed  equal 
(P(H  ) = P(Hj))  ; hence  (17)  can  be  written  as 

PNm.N,  [<*m  - X,t  ) (XPi  -X,))>  PNm<Np  UXm  - X,  ).  (Xp  — X,  ) ] . (18) 

The  probability  distributions  are  assumed  to  be  Gaussian.  Even  if  they  are  not,  this  assumption 
often  yields  good  tests.  The  test  can  then  be  written  equivalently  as 

( V,  - F ) ' R,7'  ( V,  - F ) < ( V - F ) ' Rj\  ( V , -F),  (19) 

J J ‘'m‘Jp  J J ' 1 * mllp  1 


\X„,  - X, 


f R V V P V \ 


K \xP,  -<iand^  i«.Cn>vv;i 

Rather  than  rewrite  the  hypothesis  test,  the  statistical  distance  is  defined  as 

D,  = (Y,  - V,)Rn~V  <K  ~ K>- 


The  inverse  of  the  covariance  matrix  can  be  rewritten  using  a vector  identity  as 

Q _ R -1  _ 011  012 

Q - ~[021  022  r 

where 


Q\\  = (Rn„n„  ~ Rnmr,„  Rn„np  Rnpnm  > '• 

ft  2 011  Rnmnf  R npnp  ’ 

021  = ~Rn„n„  Rn„nm  011’ 


022  = R„;l  <t  + Rnnnm  0|  1 RJ„  > 


V 
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The  statistical  distance  defined  in  (20)  is  quite  complex,  and  many  of  ihe  parameters  in  the  lest 
are  unknown  As  before,  the  unknown  parameters  in  the  test  must  be  estimated  Ihe  best 
least-square  estimate  V of  the  target's  true  position  given  the  predicted  and  measured  posi- 
tion by  standard  calculations  is 

•V,  = A,  - (Oil  + Oss  + 0,2  + (*>,  > 1 1(0,,  - 0,2  )X„  + <0:2  ~ 0:i  > 1.(22) 

The  other  quantities  to  be  estimated  are  as  follows  The  covariance  matrix  Rn  „ is  estimated 
at  tne  radar,  and  the  covariance  matrix  R„  „ is  obtained  from  the  tracking  filter  For  rather  ra- 
pid update  times  the  crosscorrelations  Rn  n and  the  mean  I can  be  estimated  by  averaging 
the  results  over  the  last  few  updates.  The  procedure  just  described  is  quite  complicated,  but  the 
statistical  distance  can  be  simplified  for  a special  case. 


The  first  simplification  is  to  set  the  mean  value  of  I to  zero,  which  means  Xm  and  Xp 
have  mean  values  equal  to  the  true  target's  position  and  are  therefore  unbiased  estimators.  The 
second  simplification  is  to  set  the  crosscorrelation  between  the  measured  and  predicted  posi- 
tions Rn  n to  zero.  The  statistical  distance  (20)  then  reduces  to 


D = (v„ 


v„>  R„ 


( A'  - A ) 


+ < Xp  ~ 


(,V„ 


A,  ). 


(23) 


w here 


X.  = IR„ 


+ R. 


-I  V -I 


(R 


n„n„  * m 


X..  + R„ 


Xp  > • 


The  quantity  D is  the  sum  of  two  squared  Mahalanobis  distances.’ 


(24) 


Equation  (23)  is  a good  statistical  distance  measure  if  the  measured  and  predicted  posi- 
tions are  unbiased  estimates,  which  is  usually  the  case  when  the  targets  are  nonmaneuvering, 
the  tracks  have  nor  previously  been  corrupted,  and  the  measurements  are  unbiased  because  of 
the  use  of  good  radar  and  radar-alignment  procedures  However,  if  this  condition  of  nonhiased 
estimates  does  not  hold,  the  mean  value  of  the  f,  and  the  crosscorrelations  between  the 
predicted  and  measured  positions  should  be  estimated,  since  they  are  now  nonzero. 


The  question  naturally  arises  as  to  how  much  better  the  test  given  by  (20)  is  over  the 
simplified  test  given  by  (23)  Is  it  worthwhile  performing  all  the  extra  computing  of  estimating 
the  mean  of  f and  the  crosscorrelation  terms.  This  question  is  not  addressed  directly  in  this 
report.  The  answer  requires  a long-term  simulation  and  analysis  procedure  with  the  use  of 
many  typical  situations.  However,  a few  comments  are  in  order  First,  given  that  the  radar 
measurements  are  unbiased  and  no  bias  exists  between  different  radars,  the  predicted  position 
is  unbiased  typically  when  the  target  is  not  maneuvering  or  the  track  has  not  been  previously 
corrupted  most  of  the  time.  The  simplified  statistical  distance  is  sufficient  under  these  cir- 
cumstances. If  the  predicted  position  is  corrupted  because  of  target  maneuvers  or  previously 
misassigned  detections,  the  estimates  required  for  the  hypothesis  test  will  be  rather  poor  any- 
way, and  proper  sorting  of  the  data  will  probably  require  a longer  time  history  of  detections. 
Therefore  it  is  felt  that  the  simpler  hypothesis  test  given  by  (23)  is  sufficient  to  use  in  all  cases. 
The  statistical  distance  measure  (23)  can  also  be  used  as  a turn  detector  in  that  if  the  detection 
and  tracks  are  properly  associated  and  the  statistical  distance  is  large,  the  target  is  turning  and 
hence  the  bandwidth  of  the  filter  should  be  opened. 


*(i  V Trunk  and  J I)  Wilson.  "Tracking  Filters  lor  Multiple-Platform  Radar  Integration.  NRI  Report  8087.  Dec  14. 
1976 
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SUMMARY 

The  concept  of  correlating  or  associating  detections  with  tracks  uses  a measure  commonly 
called  statistical  distance.  With  the  use  of  basic  concepts  in  hypothesis  testing,  the  statistical 
distance  commonly  used  in  search-radar  tracking  systems  was  derived.  By  use  of  more  of  the 
information  available  that  is,  by  use  of  predicted  positions  as  well  as  measured  positions,  a new 
hypothesis  test  was  constructed  and  a new  statistical  distance  was  obtained  which  would  prob- 
ably be  of  more  use  in  integrating  radar  data  from  noncollocated  sites.  In  general  the  form  is 
quite  complicated,  requiring  some  rather  difficult  estimation.  Under  a special  case  the  statistical 
distance  becomes  easy  to  implement  and  is  essentially  the  same  as  the  Mahalanobis  distance. 
All  that  is  essentially  required  is  that  the  measured  and  predicted  positions  be  unbiased  estima- 
tors, which  is  true  in  most  cases.  In  the  cases  this  is  not  true,  such  as  maneuvering  targets,  cor- 
rupted tracks,  and  biased  radar  data,  it  is  generally  thought  that  the  general  form  of  the  statisti- 
cal distance  is  not  neccessary  to  compute  because  other  factors  must  enter  into  the  correlation 
process  as  well.  Finally  the  statistical  distance  can  be  used  to  determine  if  the  target  is 
maneuvering. 


